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Abstract 



In |BI01 | we have proven that the generating function for self-avoiding branched 



polymers in D + 2 continuum dimensions is proportional to the pressure of the hard- 
core continuum gas at negative activity in D dimensions. This result explains why 
■ the critical behavior of branched polymers should be the same as that of the itp^ (or 

Yang-Lee edge) field theory in two fewer dimensions (as proposed by Parisi and Sourlas 
in 1981). 



In this article we review and generalize the results of [BI01]. We show that the gen- 
erating functions for several branched polymers are proportional to correlation func- 
tions of the hard-core gas. We derive Ward identities for certain branched polymer 
correlations. We give reduction formulae for multi-species branched polymers and the 
corresponding repulsive gases. Finally, we derive the massive scaling limit for the 2- 
point function of the one-dimensional hard-core gas, and thereby obtain the scaling 
form of the 2-point function for branched polymers in three dimensions. 

Keywords: branched polymers, Yang-Lee edge, repulsive-core singularity, dimensional re- 
duction, hard rods 



♦Research supported by NSF Grant DMS-9706166 



1 Introduction and Main Results 



We define the generating function for branched polymers (mod translations) to be 



00 r,N r 
N=l T JM 



D+2)(iV-l) 



dy2...dy N l[[2U , (\y i -y j \ 2 )]l[U(\y i 

ij€T ij<£T 
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Here y\ — 0, y^ . . . , are the positions of the monomers, and we sum over all tree graphs 
T on {1,... ,N}. We assume that U{t) is a positive weight function which tends to 1 
as t — > 00, and that U'(\y\ 2 ) is a positive, integrable function of y G M D+2 . By taking 
limits, we may take U(t) = 9{t — 1), where 9 is the Heaviside step function. In this case, 
2U'{\yi — Vj\ 2 ) = 5{\yi — Vj\ — 1), and we obtain our standard model of hard spheres such 
that spheres i and j are required to touch if ij G T. 

The above definition is a direct translation to the continuum of the familiar model of 
lattice branched polymers. On the lattice Z* D+2 , a branched polymer is a finite connected 
set of nearest-neighbor bonds with no cycles [[Sla99|| . An N- vertex branched polymer is a 
subset {y±, . . . , |/at} of Z D+2 , together with a tree graph on {yi, . . . , y^} such that for every 
{ViiVj} £ T, \yi — yj\ = 1. One defines cn to be the number of iV-vertex branched polymers 



mod translations. Then, as in [ Fro86 |, the generating function Zbp(z) = ^2 N z c N can be 
written as 



N 



N=l T V2,-,y N ij£T ij£T 
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where 2f//- 



Su 



. y .\ t i and t/y 



enforce the adjacency and loop- free conditions, 



respectively. For example, C3 = 6 in Z 2 , which is correctly accounted for in Q1.2 ) as there are 
3 trees, 4 possibilities for 1/2, and then 3 for y%. For more details, see |BI01| . 

Returning to the continuum, we define the partition function for the repulsive gas in a 
box A C 



°° Z N f 

Zhc(z) = / dx 1 ...dx N Yl U (\ 

N=0 ' ^ AN l<i<j<N 
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The main result of [BI01] is that the identity 



li ~D~\A\ l ° gZuC ^ 



:i.4) 



holds for all z such that the right-hand side converges absolutely. The left-hand side of 
( |1.4| ) is l/(kT) times the pressure of the repulsive gas. Evidently, its leading singularity 
~ (z — z c ) 2 ~ QHC is identical to the leading singularity ~ (z — z c ) 2 ~ 7BP of Z^p, where z c is the 
closest singularity to the origin. Hence 



ohcP) =7bp(£> + 2). 



1.5) 



2 



If one can define 9 from the asymptotic form cn ~ z~ N N~ e , then 6 = 3 — 7hc by an Abelian 
theorem. Furthermore, one expects that a, the Yang-Lee edge exponent, is equal to 1 — «hc 
Pf '99 , eqn. 7], which with ( |1.5|) leads to the Parisi-Sourlas relation | PS81 | 



6{D + 2) = a(D) + 2. (1.6) 

One can also see that the exponents u BP , r] BP are equal to their hard-core counterparts in 
two fewer dimensions (see Section |3[). 
If one takes 

U(\xi -x 3 \ 2 ) = e ~ w{x ^\ (1.7) 
with w(k) > 0, then by the sine-Gordon transformation, (|1.3|) can be written as 



^hc(z) = j exp Udxze*^ d[i w (ip), (l.{ 



where dp w is the Gaussian measure with covariance w, and z = ze w (°^ 2 . Thus certain 
branched polymer models can be written as — ze tip field theories in two fewer dimensions. 
Taking into account an effective mass term ~ <f{x) 2 from dp w , one finds a critical value for 
if on the imaginary axis, and at the critical z the interaction is itp 3 + higher order. Thus 
one expects that these theories are in the same universality class as the Yang-Lee edge (i(f 3 ) 
theory |Fis78|| (see [LF95| , |PF99|| for a more complete investigation of the hypothesis that 



the repulsive-core singularity is in the Yang-Lee class). We note that Shapir [ Sha83| , Sha85 



has given a field theory representation for lattice branched polymers which reduces to the 
supersymmetric Yang-Lee model of [|PS81|1 when presumably irrelevant terms are dropped. 



Cardy has argued recently ||Gar01|| (see also his contribution to this issue) that the 



crossover from area-weighted self-avoiding loops to ordinary self-avoiding loops in two di- 
mensions is governed by a scaling function related to the Airy function. Part of his argument 
is the reduction of two-dimensional branched polymers to zero-dimensional itp 3 theory. This 
is, in essence, the content of equations ( |1.4| ) and ( |1.8| ). 

Correlation Functions 

We define first the basic n-point density correlations for branched polymers and for 
repulsive gases. Let 

N N 

p(x) = ^2S(x -Xi), p(y) = ^5(y-yi), (1.9) 

i=l i=l 

where x, Xi G MP and y, yi G M' D+2 . Then we put 

G^kyi, ■ ■ ■ , in, z) = E m E / d yi ...dy N \{ pim) \{ [2U> 3 ] n Z7 y 

N=l ' T J^ (D+2)N i= i ij € T ij£T 

Chc(xi,... ,x n ;z) = hm D ( J[ ) ■ C 1 - 10 ) 

\i=l / HC,A 
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Here U-j := U'(\yi — yj\ 2 ), U%j : = Ui\y% — Vj\ 2 ), and (-)hc,a is the expectation in the measure 
for which Zhc( z ) is the normalizing constant. We also write G^ T for the corresponding 
truncated expectation. 

If yi, . . . ,y n are distinct points, then we have 



G^(y u . . . , y n - z) = m E / ■ ■ ■ U U U *r C 1 ' 11 ) 

M=0 ' Tan{l,...,n+M} ij&T ij<£T 

Thus, for distinct points G^p is a sum/integral over branched polymers whose vertices include 
yi, ■ ■ ■ ,y n - When points are not distinct, Gpp and G^q are understood by smearing each 
p(y) or p(x) by test functions. Thus in general, G^p and G^q are distributions which contain 
5-function singularities at coinciding points. In addition, if U[a is not smooth (for example 

in the hard sphere case U(t) = 9{t - 1)), then G$ will inherit singularities from U[y 

The density correlations G^q arise naturally when taking an order n variational derivative 
of Zhc with respect to an external field. However, we will need a different set of Green's 
functions for the repulsive gas. Stripping G^ of its singularities at coinciding points, we 
write 

^(xi, . . . ,x n ;z) = Zhc(z)" 1 E F dx n+1 . . . dx n+m U^. (1.12) 

m=0 l<i<j<n+m 

(n) 

Here n particles are forced to be at x±, . . . ,x n , and if these are distinct points, then g^c = 
C?hc- In general, G^q is a sum of terms, each with some g-gg, j < n multiplied by (^-functions 
in some of the x^s. For example, 

G n %(x 1 ,x 2 ) = g ( £ ) c (x 1 ,x 2 ) +gSc(x 1 )8(xi -x 2 ). (1.13) 

We shall see that g^ can be related to a certain n-tree branched polymer correlation func- 
tion: 



00 z p+ n r 

* * W. * ' hs(D+2) v 



Here is a loop-free graph or forest on {1, . . . ,n + p} which consists of n connected 
components or trees, each of which contains one of yi, . . . ,y n . 

Theorem 1.1. If z is in the interior of the domain of convergence at Zpp, then in the limit 
A / R D , 

g^{xi, ... ,x n ;z) = (-27r) n 5(J5 hri, . . . , x n ; , (1-15) 
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and 



G 



(n),T 
HC 



{Xl, ■ ■ ■ 



-2ir) 



dz2 ■ ■ ■ dz n Gpp | xi, 1/2, 



2tt 



;i.i6) 



)D+2 



i/ere G M. D and jji = (Xi,Zi) G 

The relation ( |1.16| ) was proven in | |B101[ | by differentiating ( |L4p with respect to sources. 
We prove ( |1.15| ) in Section |2|. As g^Q and G^q agree at non-coinciding points, ( |1.15| ) and 
( |1.16| ) combine to give relations between g^p and Gpp. In particular, we show that the 



two-point functions obey a Ward identity 

d 



d(r s 



(2) 
9bP 



-G 



(2) 
BP' 
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where by rotation and translation invariance g BP and Gp P can be thought of as functions of 
r 2 = 1 2/1 - 2/2 1 2 only. 

In an appendix, we generalize the above results to repulsive gases and branched polymers 
with more than one species of particle/monomer and species-dependent interactions. Exam- 
ples include the Widom-Rowlinson model of penetrable hard spheres |[WR70|| . As with the 
models discussed above, dimensional reduction is actually a consequence of an underlying 
supersymmetry of the branched polymer model. This requires that the attractive interaction 
between neighboring monomers be related to the repulsive interaction. 

In Section |3] we focus on the case D = 1 and derive a number of results for the standard 
hard-core gas using the method of Laplace transforms. We give fairly explicit formulas for 

(2) T ^ 1 

Gjjc and derive the values «hc = § , ^hc = |, Vac = — 1, thereby obtaining the same 
values for 7bp, i^bp, ?7bp in three dimensions. (Note that for D = 0, log^nc = log(l + z), so 
that the two-dimensional Zbp(z) has a logarithmic singularity at z = 2tt, which implies that 
7bp = «hc = 2. Unfortunately, dimensional reduction gives no information on z/ B p, t/bp in 
this case.) We derive the scaling form of the two-point function near z c = 



,(2) 



-1. 



G^ T (0,x;z) ~ \x\^ D - 2+ ^K uc (x/0, 



with 



which implies 



with 



x- 



Gg(0,y;z) ~ \x\~ d ~ 2+ ^ K BP (x/C), 



1.18) 



1.19) 



1.20) 



ir 2 x 



1.21) 



The form of Khq(x) is the same as that of the one-dimensional Ising model near the Yang-Lee 
edge [ Fis80|| . The form of Kbp{x) agrees with the prediction of Miller [ [Mil91 . 
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2 The Forest-Root Formula and Dimensional Reduc- 



tion 



We wish to derive relationships between the hard-core Green's functions in D dimensions 
and the branched polymer Green's functions in D + 2 dimensions. The key is the Forest-Root 
formula, proven in |B101|. Let f(t) be any smooth function of variables 



t 



1 < i < j < N and U 



1 < i < N, 



where each Zj £ C. Assume that /, when regarded as a function of z 1; . . . , z^v, has compact 
support. Any subset of the bonds {ij\l < % < j < N} forms a graph on the vertices 
{1, . . . , N}. A subset R of vertices is called a set of roots. A forest F is a graph that has no 
loops. The connected components of a forest are trees. We are declaring that a graph with 
no bonds and just one vertex is also a tree. See Fig. 1. 




tree 



tree 



tree 




tree 



forest 



Figure 1: Example of a forest 



Theorem 2.1. (Forest-Root Formula) 




where f( F ' R \t) denotes the derivative with respect to the variables t^ with ij £ F and ti with 
i £ R. The sum is over all forests F and all sets R of roots with the property that each tree 
in F contains exactly one root from R, and d 2 z = du dv, where z — u + iv. 

Proof of Theorem \l.l\ ( \1.13j ). In order to examine g^, we set N = n + m and put 

m n 

/(t)= n ^(Kf+yn^n^-A). ^ 

l<i<j<N i=l j=l 
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where Xij = Xi — Xj, e > 0, and where g is a smooth, decreasing function with compact 
support such that g(0) = 1. Working in a finite box A, we write 

^hc(zi, • • • , x n )Z KC (z) = ^ r / • • • dx n+m f(0), (2.3) 

m=0 m ' ™ 

and insert (|2.1p. With ?/, = (xj,Zj), this becomes 

£ ^ 1 ■ ■ ■ £ i (5) II' II' II II < 2 - 4 > 

The forest F may be divided into F^ n > (which consists of all the trees containing any of the 
vertices 1, . . . ,n), and the rest, F. If any of the vertices 1, . . . , n are not roots, then the 
corresponding factor g(t/e) is not differentiated, and so lacks a factor e~ l to compensate for 
the 0(e) integration volume. Therefore, neglecting 0(e) terms, — T\ U ■ • • U T n with 
Xi G Tj for i — 1, . . . , n , and Tj disjoint. Let m = p + M, where M is the number of vertices 
in F, and rewrite (|2.4j) as 




where o(l) denotes a quantity which tends to zero with e. Here i? is a subset of {n + p + 
1, . . . ,n + p + M} and as before, each tree of F has exactly one root from R. We have 
eliminated the integrals over z\, . . . , z n because (— 7re) _1 </(tj/e) tends to S(zi). The factors 
g(eti) with n+l<i<n + p can be replaced with 1 because the decrease of U' in essence 
forces the corresponding z^s to remain bounded. Any errors from these approximations are 
o(l). 

The only barrier to writing ( |2.5|) as a product is the presence of interactions C/y linking 
and F. However, for small e the trees in F are rarely close to each other or to the 
trees of F^ n \ Using again the decrease of U' we see that all the vertices of a tree are in 
a bounded cluster and that g(eti) can be replaced with g(et r ), where r is the root of the 
tree. Then the sum over roots leads to a factor N(T), the number of vertices in T. Observe 
that (—e/7r)N(T)g(et r ) N( - T ^~ 1 g'(et r )d 2 z r is a probability measure which becomes very wide as 
e — > 0. Thus with high probability, the interactions between F^ and F can be replaced 
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with 1, with an additional contribution to the o(l) error. As a result, ( |2.5| ) can be rewritten 

as 



P[ J (AxC^ =n " 



ff n - n * 



— 71 

p=U * " «=n+l ij^FW ij&F( n > 

Taking the limit as e — »• 0, we obtain a relation for finite .D-dimensional volume A: 



^HC 



'z) + m. 



(2.6) 



(n) / \ (n) 



(2.7) 



The limit A /* M 15 exists for each term in the sum over p, by monotone convergence. By 
dominated convergence, the sum on p may be interchanged with the infinite volume limit, 



and we obtain the first part of Theorem [LI 



Proof of ( \1-4D - A similar factorization occurs in Zhc( z ), so that all the terms with k trees 
can be written as l/k\ times the k th power of 



oo jy 
Z 



N=l ' T 



A 



N 



II II' u* 



(2.8) 



so that, as argued in [BI01 



lim. rn-log^Hc(z) 



(2.9) 



Proof of Theorem \1.1\ ( \1.16j ). The relation between (?hc and Gp^p may be derived by 
differentiating ( |2.9| ) with respect to sources. Then, as explained in [BIOl . 



\d 2 Zi G^l { 



i=2 



Z 

2^ 



(2.10) 



which is ( 1.16 ). In momentum space, then, Ghc ma y De obtained from Gpp by setting the 
components of momenta in the two extra dimensions to zero. This contrasts with the relation 



between g^ and g^p, in which the spatial components in the two extra dimensions are 
set to 0. 

Proof of ( \1.17\) . Relations between Gbp and #bp may be derived by combining (|2.7|) and 
( p.lO| ). For example, consider the 2-point functions, which by rotation invariance can be 
expressed as functions of the squared-distance t: 



» 



G@(f;z) 
4 2 p (t;z) 



Gj± c (xi, x%; z), where |xi — X2I 2 = t, 



G^l(yi,y 2 ;z), where \y t - y 2 \' 2 = t, 



(2.11) 
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and similarly for g^ c and g^p- Since g^ c agrees with G^ c at non-coinciding points, (|2.7| ) 
and (p.lO| ) imply that 

/oo 
Trdt' G^l(t\z) } t^O. (2.12) 

Differentiation yields 

j/^kt-z) = \Gtl{t-7). (2.13) 

This may be thought of as a Ward identity for the supersymmetry of our model of branched 
polymers. 

3 Green's Function for the Hard-Core Gas in One Di- 
mension 

Laplace transforms can be used to give fairly explicit formulas for the Green's function 
for one-dimensional gases with only nearest neighbor interactions. We follow [|FW69|1 in 



deriving the relevant expressions for the basic hard-core gas with no interactions other than 
a minimum separation of 1 between particles. 

Let us write the grand canonical partition function in the following way (we omit the 
subscript HC in most of this section): 

Z(L) = J2 zN dx i dx 2 --- dx N . (3.1) 

7V=0 J x±>l J X2>xi + 1 J L— l>i;jv>x']v_i+l 

The particles are restricted to the interval A = [1,L — 1], as if external particles had been 
placed at and L. We assume L > 1 and put Z(L) = 1 for 1 < L < 2. The Laplace 
transform can be evaluated explicitly: 



Z{s) := j dLe~ sL Z(L) 

oo 

= £z"J0O" +1 , (3-2) 



where 



N=0 



fOO 1 

/ dxe- sx = -e 



s 

Using analytic continuation as necessary to define Z(s), we have 



Z(s) = - J{S l , = . (3.4) 

w 1 - z J(s) se s - z v ' 
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Figure 2: Graph of the function z = se 



We obtain Z(L) by inverse transform: 

Z(L) = / -j— e sL ds. (3.5) 
2m J se s — z 

This leads to a residue formula 

n=0 

where {s n } are the solutions to se s = z, arranged in order of decreasing real part. These 
solutions are the branches of the Lambert W^-function ||CGHJK|j . 

We will make use of some properties of the s n . For z > 0, there is one real solution, and 
for — e^ 1 < z < 0, there are two real solutions (see Fig. 2). The complex solutions come in 
conjugate pairs, and all have real parts which are less than the real solutions. (This can be 
seen by writing s = x + iy and letting x(y) solve the modulus equation (x 2 + y 2 )e 2x = |z| 2 . 
Then = — l/[2(x+x 2 +y 2 )} < — c < with c independent of x, y in any bounded region not 
intersecting {(x, y)\x G ( — 1, 0)}. This shows, in fact, that the upper gap Re(si — s 2 ) > B > 
with B independent of z in any interval [—1, z ] with z < 0.) If we put s n = x n + iy n , then 
\y n — nir\ < const ||CGHJK| , Fig. 4]. In addition, x n ~ — log \y n / z \ ~ — log |?V Z | for large n, 
from the modulus equation. Hence the sum in (|3.6| ) converges for all L > 1. 

The density, or one-point function, G^ (x) = is the expectation of p{x) = J2f=i 
Xj) in the limit as A /* oo. If we take A = [— |f + 1, ^ — l] , then 

Only the n = term of ( |3.6| ) survives the L — > oo limit. Thus we have 

= = — r> (3.8) 

s + l s + l' V ; 
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where we have used the relation z = soe s °. If we identify with the density p and solve for 
the pressure p = s^kT, we obtain the equation of state for hard rods of unit length [ Ton36 1 : 

pkT 



V 



l-p 



(3.9) 



The density-density correlation, or two-point function G^ 2 ^(0,x) = G™(x), is the expec- 
tation of p(0)p(x). For x > 1, this can be written as 



G (2)( 



X) 



lim 

L— >oo 



Z (%) zZ(x)zZ - x) 



Z(L) 



so 
so + l 



-SQX 



zZ(x) 



One can insert the formula 

Z(x) = J2^(x-N-l) N 6(x-N-r 



' z N 



N=0 



(3.10) 



(3.11) 



to obtain the long- known expression for (see, for example, | SZK53| , eqn. 32]) which is 
useful if x is not too large. Alternatively, one can insert (|3.6|) to obtain 



^ J ^ S o + 1 8 n + 1 

n=0 



o (s n -s )x 



(3.12) 



and after subtracting G^ 2 , we obtain an expression for the truncated Green's function 



G (2),T, 



■T) 



s 



J s + 1 s n + 1 



^ e (s„-s )a: 



(3.13) 



which is convergent for x > 1. It is apparent from (|3.11|) that (^) 7V G^ 2 )' T (x) is continuous, 
except for a jump at x = TV + 1. This is reflected in the divergence of the series ( |3.13| ) at 
x = N + 1, when differentiated N times. For \x\ < 1, G^ 2 ^ ,T (x) = except for a 5-function 
at with coefficient G^ 1 '. 

Scaling Form of the Green's Function 

As z \ z c = — e" 1 , the two real solutions to se s = z approach the value —1. The 
correlation length is given by 



lim -— logG (2) ' T (x) 



\x — »oo \X\ 



-1 



,s ~ Si 



(3.14) 



It is clear from Fig. 2 that so — si ~ (z — z c ) a . Hence, as z \ z c , £ diverg 
with a correlation exponent z/ = |. 



es as z — z r 
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If we let z \ z c and x — > oo while keeping x = x/£ fixed, then the asymptotic form of 
G^' T is described by a scaling function 



K(x) = lim x D ~ 2+r 'G^' T \ 



x). 



x — >oo,z \z c 



(3.15) 



Here D = 1, and we take the anomalous dimension rj = —1 in order to get a nontrivial limit. 
From ( |3.13| ), and the uniform gap between si and the other solutions, we have 



K{x) 



A short calculation shows that 



lim 



S Si 



x^oo,z\z c Sq + 1 Si + 1 X^ 



so gi 



so+1 si+1 (so-«i) : 



■(1 + 0(z — z c )), and hence 



(3.16) 



AT(x) = -- e"*. 
We may also define a scaling function for branched polymers: 

K BP {x) = lim _ x d - 2+r?BP G^(0, 



x— >oo,z/z c 



Here x = x/£bp with 



£bp := 



1 



lim —f-r\ogG%l(0,x) 

»oo \x\ 



(3.17) 



(3.18) 



(3.19) 



and t]bp is chosen so as to obtain a nontrivial limit for Kp,p- As explained in [|BI01|| , we have 
the relation ( [2.10|) between G BP and Ghc > wmcn when differentiated yields 



p( 2 ) if- -—\ - —- 



G^ T (t,z). 



(3.20) 



Hence, the critical activity z c for branched polymers is equal to — 27rz c = 2-rce, £bp ( — ^) 
£hc(z),?7bp = Vkc, and 



K BP {x) 



4tt 2 
1 

7T 2 X 



[xK^ c (x) — (D — 2 + m c)Knc(x)} 



(3.21) 



Of course, the exponent z^bp governing the divergence of £bp as z f z c must equal z^hc 
if d = D + 2 = 3). 



12 



Appendix: Multispecies Examples 



We can generate the arguments of Section Q to multispecies examples. Define a repulsive 
gas partition function in D dimensions: 

oo 1 N „ 

Z HC (z) = ^ M E Il z « dx 1 ---dx N l[U ai ^(\x i -x j \ 2 ), (A.l) 

N=0 ' a lt ...,a N i=l AN ij 

where each ccj is summed over the set of species of the problem, z a is the activity of species 
a, and U alS is a repulsive interaction between species a and species j3. The corresponding 
multispecies branched polymer generating function is 



A? 



^ i r 

^ = e e e n z « L +2W n to n ^ ( a - 2 ) 

where 

Ui^U^dxi-Xj] 2 ), (A.3) 

and U-j is its derivative. In particular, for the hard-core model with minimum separation 
R a/ 3 between species a and species /3, we would have U al3 (t) = 6{t — R^g) and 

2U al3 \\x i - Xj \ 2 ) = -^-5(\ Xi - Xj \- R aP ). (A.4) 



ai3 



Assume that U a ^ satisfies the usual conditions (U af3 , U af3> positive, U al3 1 at oo, U a/3 ' 
integrable in M £)+2 ). Then, provided Z B p(z) is absolutely convergent, we obtain a reduction 
formula 

lim JL log Zhc(z) = -2nZ BP . (A.5) 



We also obtain results as in Theorem |1.1| for correlation functions. 
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